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Abstract 


We searched for superflares on solar-type stars using Kepler data with 1 min sampling in order to detect 
superflares with short duration. We found 187 superflares on 23 solar-type stars whose bolometric 
energy ranges from the order of 10^^ erg to 10^® erg. Some superflares show multiple peaks with the 
peak separation of the order of 100-1000 seconds which is comparable to the periods of quasi-periodic 
pulsations in solar and stellar flares. Using these new data combined with the results from the data 
with 30 min sampling, we found the occurrence frequency (dN/dE) of superflares as a function of flare 
energy {E) shows the power-law distribution {dN/dE oc with a ^ —1.5 for 10^^ < E < 10^® erg 
which is consistent with the previous results. The average occurrence rate of superflares with the energy 
of 10^^ erg which is equivalent to XlOO solar flares is about once in 500-600 years. The upper limit of 
energy released by superflares is basically comparable to a fraction of the magnetic energy stored near 
starspots which is estimated from the photometry. We also found that the duration of superflares (r) 
increases with the flare energy {E) as r oc i50-39±o.03_ explained if we assume the time-scale 

of flares is determined by the Alfven time. 
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Introduction 

Solar flares are eruptiv e events in the solar atmosphere caused by the magnetic reconnection (e.g. 


Shibata &: Maeara 


2011 b The occurrence frequency of solar-flares decreases as the flare energy increases. 


The frequency-energy distribution of solar flares can be fitted by a simple power-la w function with an 


index of —1.5 - —1.9 in the flare energy range between 10^"^ erg and 10^^ erg (e.g. 

Aschwanden et al. 

2000 

Shimizu 

1995 

Crosbv et al. 

1993 

). The total bolometric energy released by the largest solar flares 


is estimated to be the order of 10^^ erg (lEmslie et al 


2012h and the occurrence frequency of such flare is 


about once in 10 years. On the other hand, muc h larger flares called ” superflares” have been observed o n 


a wide variety of stars including solar-type stars (ILandini et al 


1986 


Schaefer 


1989; 


Schaefer et al 


2000 ). 


Recently, many superflares on solar-type stars (G-type main sequence stars) have been discovered by the 
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Kepler space telescope (jMaeliara et al 


2012 ; 


Shibavama et al 


2013; 


Wu et al 


2015 ). The bolometric 


energy released by the superflare ranges from 10^^ to 10^® erg which is 10-10^ times larger than that of 
the largest solar flares (10®^ erg). Although most of G-type main sequence stars exhibiting superflares are 
rapidly rotating stars, the temperature and rotation period of some superflare stars are close to those of 


the Sun. For example, 


No garni et al. 


(|2014ll reported that the rotational velocities of two superflare stars, 


KIC 9766237 and KIC 9944137, are comparable to that of our Sun and no hint of binarity was detected. 
The average strength of the magnetic field on the surface of these stars, by using the absorption line of 
Ca 11 at A 8542 A, are estimated to be 1-20 G whi ch is comparable to or one order of magnitude higher 


than that of the Sun 2 G). On the other hand, (IWichmann et al 


201411 reported that some superflare 


stars with photometric period of 10-20 days have large rotational velocities of vsini ^ 10-20 km/s. This 
suggests that these stars are either sub-giant stars or binary systems, not main sequence stars like our 
Sun. 

Solar flares are sometimes accompanied by coronal mass ejections (GMEs). Earth-directed CMEs can 
cause geomagnetic storms which can produce geomagnetically induced currents (GlGs) in long electrical 
conductors (e.g. power lines) and can damage electrical equipment. For example, on March 13, 1989 about 
6 milli on people in Queb ec, Canada suffered an electrical power blackout caused by a huge geomagnetic 


storm (Allen et al 


l|Carrington 


3 ) 


19891) . Moreover, in 1859, one of the largest solar flares, call ed the Carrington flare 


185911 . caused the largest geomagnetic storm in the past 200 


years (Tsurutani et al 


2003) 


186111 . Even the solar 


which led to a failure of telegraph systems in Europe and North America (ILoomisI l 
flares which release up to 10^^ erg of energy can damage not only our power grids, but also communication 
and navigation systems. If superflares were to occur on the Sun, they would cause much more serious 
damage to our civilization. Therefore it is important to study the properties of superflares on solar-type 
stars. 

Previous studies of superflares on the basis of the Kepler data mainly focus on superflares with extremely 
large energy (> 10^"^ erg) and long duration (> 1 hour) and could not discuss the statistical properties 
of superflares with energy between 10®® erg - 10®'^ erg due to low time-resolution of the data. In this 
paper, we report the results from the data with high time-resolution 1 min) and discuss the statistical 
properties (e.g. occurrence frequency, flare duration) of superflares on solar-type stars. 
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Methods and Data 


We searched for flares (sudden brightenings of the star) from the short-cadence data 1 min interval) 


Koch et al. 

2010 ; 

Gilliland et al. 

2010 ) 


2010l) . We selected solar-type (G-type main sequence) stars from the 


Huber et al. 

(2014 

) instead of those from 

Brown et al. 

papers ( 

Maehara et al. 

2012 

Shibavama et al. 

2013 

), ' 


(l2011ri (initial Kepler Input Catalog) . In previous 


2013 ), we used stellar parameters taken from 


Brown et al 


(12011 ) and the selecti on criteria of 5100K < Tgff < 6000 K and log 5 > 4.0. However the te mperatures 


Brown et al 


(|201ll) are systematically lower by ~ 200 K than those in 


Huber et al 


occurrence frequency of superflares depends on the surface temperature (e.g. 


(2014 ). Since th e 


Candelaresi et al. 


20141) 


this systematic difference in surface temperature of stars would affect the resultant flare occurrence 
frequency. Therefore, in order to reduce the difference in the occurrence frequency of superflares caused 
by the systematic difference in temperature, we used the selection criteria of 5300K < Teff < 6300K and 
4.0 < logg < 4.8. The total number of solar-type stars observ ed with short-cadence mode is 1547. 


We used the Presearch Data Conditioned (PDC) light curve (IStumpe et al. 


2012 : 


Smith et al.l 


2012 l) for 


the detection of flares. The typical length of each continuous observation in short-cadence mode is about 
30 days. We first calculated the distribution of the following parameter (AT^) for all consecutive data 
points from two pairs of the stellar brightness {F(t)) at the time of ti, and (i and n are 

integer numbers and ti means time of the i-th data point): 


AF^(ti^n) = s(Fi - Fi_„_i)(T;+i - Fi_„), 


( 1 ) 


where 


s = 


( 2 ) 


1 (Fi - > 0 and F^+i - Fi^n > 0 ) 

-1 (Fi- Fi_n-i < 0 or Fi+i - Fi^ri < 0 ) 

In the case of typical 30-day short-cadence data blocks, i ranges from 1 to ~ 44000. If the flare peak 

is around at ti and the stellar brightness at is close to the quiescent level, the value of AF^{ti) 

becomes much larger than the mode of the distributi on. Therefore tj — should be comparable to or 

longer than the rise time of superflares. According to 


Shibavama et al 


on KIC 11610797 is about 2 min. In the case of the X17 [(polar flare (|Kopp et al 


I 


201, si 


( 2013D. the rise time of superflare 


2005 ). the rise time is 


about 4 min. Thus we used the separation between two data point n = 3 {U — ti-n ^ 3 min) and n = 10 
(^ 10 min) for the analysis. In the case of larger n value, we can detect flares having longer rise-time. 
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The threshold of the flare was determined to be three times the value at the top 1% of the distribution. 
We removed long-term brightness variations around the flare peak by fitting a quadratic function. We 
used the data from 0.05 to 0.01 (in the case of n = 3) or from 0.15 to 0.03 (n = 10) days before the flare 
peak and the data from 0.05 to 0.10 day (n = 3) or 0.15 to 0.25 day (n = 10) after the peak for the 
fitting. After removing the long-term light variations, the flare start time and flare end time are defined 
to be the time at which the residuals of stellar brightness become larger and smaller than three times of 
the photometric error (dphoto; PDC flux errors given in the Kepler data) respectively. We analyzed only 
flares with the duration longer than 5 min (at least 5 consecutive data points exceed Suphoto) since it is 
difficult to distinguish between true flares and false events in the case of a single data point above the 
threshold (e.g. spikes in figure 1 (h)). 

We excluded the flare candidates which have longer rise time than decay time. Since the number of stars 
observed in short-cadence mode is only ~ 1500, the possibility that two or more stars exhibit flares at the 
same time is negligible. Thus we also removed the pairs of flare events on different stars which have the 
same flare peak time. After removing these false events, we checked all light curves of flare candidates 
by eye and removed non-flare events, since our flare detection method misidentified some eclipses or 
light curve jumps as flares. We also examined the pixel level data of flare stars in order to exclude the 
contamination of flares from nearby stars or other false events. If the spatial pixel-distributions of the 
brightness of a star on CCDs are different between during the flare and at the quiescence, t he candidate of 
a flar e is revealed to be a light variation of nearby source or a non-astrophysical event (e.g. 


Maehara et al 


201211 . We removed 5 events in this step. 


We checked the Kepler Eclipsing Binary Catalog (IPrsa et al 


2011 


Slawson et al. 


2011 


Kirk et al 


2013) 


and removed all flare candidates which occurred on the eclipsing binary (14 stars). We also removed al l 


flare candidates on the stars which are listed as visual or spectroscopic binaries in 


Notsu et al 


( 2ni5a l 


(2 stars). In total, we eliminated 152 flare candidates on these 16 stars. Since the fraction of known 
binary systems wou ld be small, there may be more unknown binary systems in the sample. According to 


Notsu et al. 


(|2015all . the lower limit of binary fraction among superflare stars based on high-dispersion 


^Solar flares are classified as A, B, C, M, or X, according to the peak X-ray (1-8 A) flux as measured by the GOES 
satellite. Peak fluxes of A, B, C, M and X class solar flares are < 10“^, 10“^-10~®, 10“®-10“^ and > 10“^ 

W/m respectively. Each X-ray class is divided into a linear scale from 1 to 9 (e.g. B1 to B9, Cl to C9, and Ml to M9). An 
X17 class flare is 17 times more energetic than an XI class flare. 
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spec troscopy is about 30 %, wh ich is not so different from the binary fraction of nearby solar-type stars 


(e-g- 


Duauennov &: Mayor 


1991 b 


Bolometric energy emitted by each flar e was estimated from th e stellar luminosity, flare amplitude and 
flare duration with the same manner as 


Shibavama et al. 


(|2013l) . The lu minosity of each fla re star were 


estimated from the effective temperature and the stellar radius taken from 


Huber et al 


i n bol ometric energy was estimated from errors in the temperature and stellar radius from 


( 2014). The error 


Huber et al 


(|2014l) . Since the error in the temperature is crT^^^r ~ 200 K (o'Tatar/’^star ~ 3-4%) and that in the stellar 


radius is cr/jatar/-^star ~ 10 %, the typical error in bolometric energy would be (T_Efia„/£’flare 25 %. 

We estimated periods of long-term light variations from light curves of each star by using the discrete 
Fourier transform (DFT) method. For period analysis, we used the long-cadence data (time resolution 
of 30 min) between 2009 September (Q2) and 2013 April (Q16). After removing linear trends in each 
quarter and the correcting the difference in mean flux between different quarters, we calculated the power 
spectrum and chose a peak corresponding to the maximum power among peaks which are 20 times larger 
than the noise level. 


Results and discussion 

We detected 187 flares on 23 solar-type stars from the data of 1547 solar-type stars. The details of all 
detected flares and properties of superflare stars are summarized in table 1-6 in Additional file 1 (CSV 
version of table 1 and table 2-6 are also available as Additional file 2 and 3 respectively). Figure 1 
shows typical light curves of superflares from short-cadence data. Light curves of all detected flares are 
shown in figure 1-94 in Additional file 1. The amplitudes of flares in a unit of average brightness of 
the star range from 1.3 x 10“^ to 8.5 x 10“^ and the bolometric energies of flares range from 2 x 10^^ 
to 8 X 10^® erg. The lower end of the amplitude of detected flare is compara ble to that of the X17 


class solar flare (normalized amplitude ^ 3 x 10 ^) occurred on Oct 28, 2003 (|Kodd et al 


200511 ■ As 


shown in figure 1 (b) and (d), some flares show multiple peaks with the peak separation of 100 - 1000 


the Sun (e.g. 


2005 


Asc 

iwanden et al. 

1999; 

Nakariakov & Melnikov 

2009 

) and M dwarfs ( 

Mitra-Kraev et al. 


Welsh et al 


200 


Anflnogentov et al 


201311 . On the analogy of the solar flare oscillations, these 


multiple peaks may be caused by coronal loop oscillations associated with flares. 
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Occurrence frequency distribution 

Figure 2 represents the occurrence frequency distribution of superflares with energy > 10^^ erg as a 
function of the rotation period. The flare frequency strongly depends on the rotation period and the 
frequency rapidly decreases as the rotation period increases in the period above a few days. No flare was 
detected on the stars with Prot > 20 day from short-cadence data and the upper limit of the frequency for 
this bin is < 1 in 80 years. The flare occurrence frequency on the stars with 2.5 < Prot < 20 day roughly 
proportional to Pj-ot^. The flare frequency on the stars with P^ot > 20 day expected from this decrease 
trend is ^ 0.004 year“^ star which is about 8 times lower than that on the stars with 10 < Prot < 20 
day. The similar relation between flare frequency and rotation period was reported by previous studies 


using the long-caden ce data (e.g. 


Notsu et al 


the age of stars (e.g. 


Ayres 


2013bll . The rotation period is thought to be related to 


19971). This result suggests that the young solar-type stars would exhibit 


superflares more frequently. 

Figure 3 (a) represent the occurrence frequency of superflares as a function of the bolometric energy 
of superflares. Solid histogram in figure 3 (a) represents the frequency distribution of superflares on all 
solar- type stars derived from short-cadence data and dashed histogram represents that from long-cadence 


data (|Shibavama et al. 


201311 . The average flare frequencies in a given energy bin from short-cadence data 


are about 5 times higher than those from long-cadence data. This difference in the flare frequency is 
mainly caused by that in the rotation period distribution of the observed targets. The period distribution 
of the targets in short-cadence mode is biased to the shorter-period end. The fraction of stars with the 
rotation period (Prot) < 10 days observed in long-cadence mode is 8.1 % (8312/102709). On the other 
hands, that of stars with Prot < 10 days in short-cadence mode is 32% (499/1547). As mentioned above, 
the stars with shorter rotation period exhibit more frequent superflares. Therefore, the average flare 
frequency derived from short-cadence data is higher than that from long-cadence data. Solid histogram 
in figure 3 (b) represents the corrected frequency distribution of superflares from short-cadence data 
taking into account the difference in the period distribution of the observed samples. Since the ratio of 
total observation time (sum of the length of observation time for each star) for stars with Prot < 10 day 
to that for stars with Prot > 10 day in short-cadence mode (1:1.7) is about 6.5 times larger than the ratio 
in long-cadence mode (1:11), we calculated the corrected occurrence frequency by adjusting the ratio of 
total observation time in short-cadence mode to that in long-cadence mode. Both frequency distributions 
from short- and long-cadence data are almost the same for the flare energy between 10^"^ erg and 10^® 
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erg. The frequency of flares derived from the long-cadence data is less than that from the short-cadence 
data for the flare energy below 10^^ erg. This difference is caused by the detection limit of flare search 


method. Since the time-resolution of long-cadence data is much lower than that of short-cadence data, 
we could not detect smaller flares from long-cadence data and observed flare frequency from long-cadence 
data becomes much lower than that from short-cadence data. 

The frequency distribution of superflares from the short-cadence data can be htted by a power-law 
function with the index of —1.4 ± 0.2 for the energy between 4 x 10^^ erg and 1 x 10^® erg. This value 
is comparable to the power-law index for the frequency distribution of superflares derived from the long- 
cadence data. Using the combined data set from both short- and long-cadence data, the power-law index 
is —1.5 ± 0.1 for the flare energy of 4 x 10®® -lx 10®® erg. 


Shibata et al 


(|201,‘1^ pointed out that the frequency distribution of superflares on Sun-like stars (early- 


G dwarfs with Pj-ot > 10 day) and those of solar flares are roughly on the same power-law line. 
Figure 4 represents the comparison between the frequency distribution of superflares on Sun-like stars 
(5800 < Teff < 6300 and Prot > 10 day) derived from both short- (Hlled circles) and long-cadence (solid 
line) data and those of solar flar es (dashed line s ). Th e thin dotted line indicates the power-law function 


with an index of -1.8 taken from 


Shibata et al 


(|2013ll . The frequency distribution of superflares derived 


from short-cadence data is also on the same power-law line. This distribution suggests that average 
occurrence frequency of superflares on Sun-like stars with bolometric energy of 10®® erg, 10®^ erg, and 
10®® erg are once in ^ 70 years, ^ 500 years and ^ 4000 years respectively. As mentioned in previous 
paragraph, the frequency of superflares strongly depends on the rotation period. The flare frequency on 
the stars with the same rotation period as the Sun (Prot ^ 25 day) would be ~ 8 times smaller that that 
on the stars with 10 < Prot < 20 day. Therefore the occurrence frequency of superflares with energy 
> 10®® erg would be once in 50 0-600 years. These v al ues are lower than t he upper limits derived from 


radionuclides in lunar samples (|Schriiver et al 


2012 ). 


Mi^ake et al 


the content of tree rings around A.D. 774-775 and 


Miyake et al 


(12012 ) reported a rapid increase in 


(2013|) found another increase 


event around A.D. 993-994. The cosmic-ray flux of these events are equivalent to solar flares with the 
energy of ~ 10®^ ergs, if cosmic-ray sources of these events were solar flares. The occurrence frequency 
of these events (two events in 270 years) is roughly comparable to the average occurrence frequency of 
superflares on Sun-like stars with the en ergy of 10®® - 10®^ erg. In the case of solar flares, the fraction of 


geo-effective CMEs is less than 10% (e.g. 


Goualswamv et al 


200711 . If the fraction is same for superflares 






























on our Sun, the occurrence rate of extreme GIC events (10 times or more larger than the Carrington 
event) would be less than once in a few thousand years. 

Relation between flare energy and area of starspots 


Most of superflare stars show large-amplitude light variations with periods from a few days to ^ 20 days. 
If we assume these quasi-periodic variations are caused by the rotation of the star , the amplitude an d 


period of light variations correspond to the area of starspots and the rotation period (|Notsu et al. 


201 3bh . 


The large-amplitude light variations of superflare stars suggest that the surface of superflare stars are 
corved by large starspots. 

Figure 5 shows the scatter plot of flare energy as a function of the area of starspot group (Agpot). The 
area of starspots on superflare stars were estimated from the amplitude of light variations (AFi-ot) by 
using the following equation: 


^spot - AFl-ot^S 


1 - 




spot 




-1 


( 3 ) 


where Agtar is the appa rent area of the st ar, Tgpot and Tgtar are the temperature of starspot and 


photosphere of the star (|Notsu et al 


2013bll . Here we define the amplitude as the brightness range 


normalized by the average stellar brightness, in which the lower 99% of the distribution of the brightness 
difference from average, except for the flares, are included. We estimated the amplitude from the long- 
cadence data of the quarter at which each superfiar es detected from th e short-cadence data occurred. 
Hst.ar and Tstar for each superflare star were taken from|j 


(|Berdvugina 


Huber et al 


(l2014ri and we assumed Tspot = 4000K 


2005(). Although there is large scatter, the majority of energetic flares occur on the stars 


with large starspots. Since the flares are a result of a sudden conv ersion of magnetic energy to other 
forms (i.e. thermal and kinetic energy; e.g. IShibata fc Magaral 1201 111 , the flare energy is supplied from 


the magnetic energy stored near the star spots (Amag) and the bolometric energy released by the flar e 


may be a fraction of the magnetic energy (IMaehara et al 


2012 


Shibata et al. 


2013 


Notsuetal 


2013bll . 


Here the stored magnetic energy can be roughly estimated by 




Stt ’ 


( 4 ) 


where B and L correspond to the magnetic field strength of the starspots and the scale length of the 
starspot group respectively. By using the total area of starspot group Agpot ~ the energy released by 
a flare can be written as 

/•n2 ^ 3/2 

Aflare /F^mag ^ (5) 


Stt 
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where / is the fraction of magnetic energy released by the flare. This suggests that the flare en- 
ergy is roughly proporti onal to the area of the starspot group to the power of 3/2. Since / = 0.1 


(Aschwanden et al 


201411 . the typical sunspot area for generating XlO-class flares 1 0^^ erg) observed 
in 1989-1997 was 3 x 10“"^ of the half area of the solar surfa ce (ISammis et al.l 1200011 and the typical 


magnetic field strength of sunspot is the order of 1000 G (e.g. 
written as 

i^flare 7 X 10^^ (erg) 


Antia et al. 


2003^ . equation (5) can be 




0.1 J VlOOOG 




A,pot/27ri?|l'/' 


0.001 


( 6 ) 


Solid and dashed lines in the figure 5 represent equation (6) for / = 0.1, B = 1000 and 3000 G respectively. 
Majority of superflares detected from short-cadence data (filled-squares) and almost all solar flares (small 
dots) are below these lines. This suggests that the upper limit of the energy released by the flare 
is basically determined by the magnetic energy estimated from the area of starspots. However the 
bolometric energies of superflares detected from long-cadence data (small-crosses) are biased to larger- 
energy end. As mentioned above, due to the low time-resolution, only superflares with larger bolometric 
energies could be detected from long-cadence data. In addition, superflares with large energy but small 
occurrence frequency could be detected from the long-cadence data because of the large number of stars 
(~ 105 ). 

There are some superflares above the analytic relation. KIC 7093428 is one of superflare stars exhibiting 
superflares whose energy is much larger than that expected from equa tion (6). This obj ect is located 


201411 . The X-ray 


within the error circle of an X-ray source ISXPS J184925.7-1-423901 (lEvans et al. 
luminosity of the object is ~ 10^® erg/s if the distance of the obj ect is ^ 200 pc. The empirical relation 


between the X-ray luminosity and the magnetic flux of the star (jPevtsov et al. 


200311 suggests that the 


magnetic flux of KIC 7093428 is the order of 10^^ Mx. This value corresponds to the spot area of 3 x 10^° 
cm^ for B = 3000G, and the energy of the largest superflare on KIC 7093428 is comparable to the flare 
energy estimated from equation (6). The spot area of the object expected from the X-ray luminosity 
is more than one order of magnitude larger than that estimated from the amplitude of quasi-periodic 
modulations. The amplitude of light variations du e to the rotation of star is affected by the inclination 
of rotation axis and the latitude of starspots (e.g. 


Notsu et al. 


2013al lbfl. In the case of stars with low- 


inclination angle or stars with starspots around polar region, the light variation caused by the rotation 
become small and the are a of starspot s derive d from the amplitude of light variation is smaller than the 


actual area. According to 


Notsu et al 


(|2015bfl , some superflare stars which show superflares with energy 
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larger than that expected from the amplitude of light variation have low-inc lination angles. 


In addition, polar spots are often observed on rapidly-rotating cool stars (e.g. 


Strassmeiei 


possibility is flares on the companion star. 


Kitzeetal 


2009 ). Another 


(|2014h found that the astrometric position of the 


Sun-like superflare star, KIC 7133671 which show a superflare with energy larger than expected from 
the photometric spot area, was shifted by 25 mas during the flare. This suggests that the flare on KIC 
7133671 occurred on a faint companion star. In this case, the flare energy would not be related to the 
area of starspot group on the primary star. 

Correlation between flare duration and flare energy 
Figure 6 shows the duration of superflares as a function of the bolometric energy of flares. Here we 
defined the flare duration as the e-folding decay time of flare intensity after its peak. The duration of 
superflares increases as the flare energy increases. The linear fit for the superflares from short-cadence 
data in the log-log plot yields 

_ _ 77>0.39i0.03 fn\ 

Tflare OC (7) 

where Tflare and Aflare indicate the duration and bolometric energy of flares. The duration of superflares 
derived from long-cadence data (crosses) is longer than that from short-cadence data. This difference may 
be caused by the difference in time resolution of the data and selection bias. Since the time resolution of 


method for long-cadence data used in 


Maehara et al. 

(2012 

) and 

Shibavama et al. 

(2013) 


superflares with long duration. The bolometric energy and duratio n in optical wavele ngth of the X17solar 


flare on 2003 October 28 are the order of 10^^ erg and ^ 10 min (jKoDP et al 


20051) . which are roughly 


comparable to the correlation between the duration and energy of superflares on solar-type stars. 
Similar correlation between the flare duration and energy was observed in solar flares. The power-law 
slope for the correlation between the duration of solar flares and X-ray fluence observed with the GOES 


is about 1/3 (|Veronig et al 


20021) . The correlation between the durat ion and peak flux observed with 


RESSI also shows similar power-law slope of ^ 0.2 (|Christe et al 


20081) . These similarity between solar 


flares and superflares on solar-type stars suggests that solar flares and superflares are caused by the same 
physical process (i.e. reconnection). 

As discussed in the previous subsection, the flare energy is related to the magnetic energy stored near 
the starspots as follows: 

Eflare ~ /E„,ag OC f , (8) 
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where /, B, and L correspond to the fraction of energy released by a flare, the magnetic field strength of 
the starspots and the scale length of the starspot group respectively. On the other hand, the duration of 
flares, especially the duration of impulsive phase of flares, is comparable to the reconnection time (Tree) 
which can be written as 

Tflare ^ Tree ^ '^AI^A ^ BjVp^j (9) 


where ta = Ljvx is the Alfven time, ta is the Alfven velocity, and Ma is the non-dimensional 
reco nnection rate which is c omparable to 0.1-0.01 for the fast reconnection such as Petschek model 
(e-g- 


Shibata fc Magara 


201lh . Since all G-type main sequence stars have similar stellar properties (e.g. 


Teff and logg), B and ta might not be so different among them. Therefore, from equation (8) and (9), 
the duration of flares can be written as 

Tflare « (10) 


This suggests that the power-law slope for the correlation between the flare duration and flare energy is 
about 1/3 and this is comparable to the observed value of 0.39 ± 0.03. 


Conclusions 

We found 187 superflares on 23 solar-type stars from the 1-min cadence data obtained with the Kepler 


(1) The power-law frequency distribution of superflares (e.g. 


Maehara et al. 

2012 

Shibavama et al. 

2013) 


continues to the flare energy of 10^^ erg and the average occurrence rate of the XlOO superflares on 
the stars with similar rotation period to the Sun is about once in 500-600 years. 

(2) Some of superflares show multiple peaks during a flare with the peak separation of 100-1000 seconds. 
The time scale of these modulation during the flare is comparable to the periods of quasi-periodic 
pulsations in solar and stellar flares. 

(3) The flare duration increases with the flare energy as Tflare oc 7?flare^°'°^- 

Observations of solar flares, especially white light flares, with the high temporal resolution are also 
important to compare the properties of solar flares and those of superflares on solar-type stars. 
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BJD - 2400000 Day from flare peak 

Figure 1. Light curves of typical superflares, (a) Long-term light variations of the G-dwarf KIC 
11610797. The vertical axis means the relative difference between observed brightness of the star and 
the average brightnes during the observation period. The horizontal axis means the times of the 
observations in Barycentric Julian Date, (b) Enlarged light curve of a superffare observed around 
Barycentric Julian Date (BJD) 2454981.63 (indicated by the down arrow in panel (a)). Fihed-squares 


18 




































with solid-lines and dashed-lines represent the light curves from short- and long-cadence data 


respectively. The long-cadence light curve is shifted in vertical direction for clarity. 

The flare peak (to) is BJD 2454981.6315 and the bolometric energy of this flare (ifflare) is 1.2 x 10^® erg. 
(c) Same as panel (a), but for KIC 11551430. (d) Same as panel (b), but for the superflare observed on 
KIC 11551430 around BJD 2456208.35 (down arrow in panel (c)). (e) Same as panel (a), but for KIC 
10338279. (f) Same as panel (b), but for the superflare observed on KIC 4543412 around BJD 
2455208.68 (down arrow in panel (e)). (g) Same as panel (a), but for KIC 11253827. (h) Same as panel 
(b), but for the superflare observed on KIC II253827 around BJD 2456365.27 (down arrow in panel (g)). 
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Figure 2. Occurrence frequency distribution of superflares as a function of the rotation 
period. The horizontal axis indicates the rotation period of the star derived from photometric data. 
The vertical axis indicates the number of superflares per star and per year in each period bin. Period 
bins are defined as follows: P < 2.5 day , 2.5 < P < 5 day, 5 < P < 10 day, 10 < P < 20 day, and 
P > 20 day. Solid-line histogram represents the occurrence frequency of superflares with an energy of 
> 10^^ erg in each period bin. The dashed down-arrow indicates the upper limit of the flare frequency. 
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Figure 3. Occurrence frequency distribution of superflares as a function of flare energy. 

(a) Bold solid and dashed histograms represent the occurrence frequency of superfiares on all G-typ e 


main sequence stars detected from short- (this work) and long-cadence data (|Shibavama et al 


2013)as a 


function of total energy of superfiares respectively. The vertical axis indicates the number of superfiares 
per star, per year, and per unit energy. Error bars represent l-cr uncertainty in the frequency estimated 
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from the square-root of the number of detected flares in each energy bin and the uncertainty in flare 
energy, (b) Same as (a), but for the corrected frequency distribution from the short-cadence data. In 
long-cadence mode, the ratio of total observation time (number of stars x observation time of each star) 
for stars with Pj-ot < 10 day to that for stars with P^ot > 10 day is about 1:11. On the other hand, the 
ratio in short-cadence mode is about 1:1.7. The corrected occurrence frequency was calculated so that 
the ratio of total observation time in short-cadence mode is same as that in long-cadence mode. 
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Figure 4. Comparison between occurrence frequency superflares on Sun-like stars and 
those of solar flares. Filled-circles indicates the occurrence frequency distribution of superflares on 
Sun-like stars (G-type main sequence stars with Pj-ot > 10 days and 5800 < Tgg < 6300K) derived from 
short-cadence data. Horizontal error bars represent the range of each energy bin. The definition of the 
vertical error bars is the same as figure 3. Bold-solid li ne represents the power- law frequency 


distribution of superflares on Sun-like stars taken from 


Shibavama et al 


the po wer-law freque ncy distribut ion of solar flares observ ed in hard X-ray (ICrosbv et al 


X-ray ([Shimizu 


I 995 I) . and EUV (jAschwanden et al 


(I 2 OI 3 ). Dashed lines i ndicate 


199.3h . soft 


2000 1. Occurrence frequency distributions of 


superflares on Sun-like stars and solar flares are roughly on the same power-law line with an index of 
— 1.8 (thin-solid line) for the wide energy range between 10^"^ erg and 10^® erg. 
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Figure 5. Scatter plot of flare energy as a function of spot area. The lower and upper 
horizontal axis indicate the area of starspot group in the unit of the area of solar hemisphere and the 
magnetic flux for B = 3000 G. The vertical axis represents the bolometric energy released by each flare. 
Filled-squares and small-crosses indicate superfiares on G-type m ain sequence stars detected from 


short-(this work) and long-cadence data (|Shibavama et al 


2 OI 3 II respectively. Small filled-circles 


represent solar flares (Ishii et al., private communication, based on the data retrieved from the website 
of the National Geophysical Data Genter of the National Oceanic and Atmospheric Administration 
(NOAA/NGDC), Solar-Terrestrial Physics Division at http://www.ngdc.noaa.gov/stp/). The area of 
starspot group for G-type main sequence stars was estimated from the full amplitude of light variations 
during each observation period in which the flare occurred using equation (3). We assumed that 
bolometric energies of B, C, M, X, and XIO class solar flares are 10^®, 10^®, 10®°, 10®®, and 10®^ erg 
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Figure 6. Scatter plot of the duration of superflares as a function of the bolometric 
energy. Filled-squares and small-crosses indicate superflares on G-type main sequence stars detected 
from short- and long-cadence data respectively. We used e-folding decay time as the flare duration. 
Dotted line indicates the linear regression for the data of superflares from short-cadence data. The 
power-law slope of the line is 0.39 ± 0.03. 
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